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ABSTRACT 

Experimental  data  for  both  the  charging  and  the  blowdown  process 
in  a  single  air  receiver  is  presented  and  interpreted  in  detail.   The 
data  indicate  that  heat  transfer  effects  cause  a  radical  departure 
from  adiabatic  behavior,  and  that  such  variances  can  be  explained 
qualitatively  on  the  basis  of  available  simplified  expressions  for  the 
state  of  a  gas  in  a  receiver.   It  is  shown  that  the  behavior  of  systems 
fulfilling  the  limiting  conditions  of  these  solutions  can  be  adequately 
predicted  and  useful  design  results  obtained.  Methods  for  determining 
heat  transfer  convective  conductances  for  use  in  the  simplified  solu- 
tions are  discussed  and  evaluated. 

The  experimental  work  was  performed  from  January  1959  through 
April  1959  at  the  United  States  Naval  Postgraduate  School,  Monterey, 
California 
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1.  Introduction 

During  blowdown  and  charging  of  a  gas  receiver,  changes  in  the 
gas  temperature  may  provide  a  substantial  temperature  difference  for 
heat  transfer  between  the  gas  and  the  receiver  walls.  In  addition 
to  heat  transfer  through  the  walls,  the  receiver  mass  itself  may  provide 
significant  energy  storage  so  that  heat  transfer  from  such  a  thermal 
mass  becomes  important.  Since  the  thermal  capacitance  of  the  receiver 
is  usually  quite  large  relative  to  that  of  the  contained  gas,  considerable 
heat  transfer  between  the  walls  and  the  gas  can  occur  with  little 
change  in  the  wall  temperature.  The  usual  approach  in  the  design  of 
charging  and  blowdown  systems  has  been  to  assume  an  adiabatic  process. 
Inasmuch  as  this  procedure  neglects  heat  transfer  effects,  substantial 
error  in  the  prediction  of  system  behavior  may  result. 

When  heat  transfer  rates  are  to  be  considered,  it  becomes  necessary 
to  specify  the  rate  at  which  mass  enters  or  leaves  the  receiver.  Two 
types  of  flow  of  considerable  interest  are  constant  mass  flow  and  flow 
through  a  critical  flow  nozzle.  The  former  case  finds  considerable 
application  in  the  blowdown  wind  tunnel.   Here  the  problem  consists  of 
maintaining  sufficiently  constant  stagnation  conditions  despite  decreas- 
ing temperature  from  the  expansion  process  in  the  reservoir.  Methods 
applicable  to  this  problem  have  been  treated  both  analytically  and 
experimentally  by  Murphy,  et  al.,  [lJ   employing  a  separate  thermal 
mass  and  air  storage  reservoir.   Incorporated  in  these  solutions  are 
analytical  expressions  for  specific  modes  of  convection  heat  transfer 
between  the  air,  reservoir  walls,  and  thermal  mass.  Thus  the  value 
of  these  solutions  is  limited  to  systems  behaving  in  a  similar  manner. 
Tfaabers  In  brackets  refer  to  references  listed  In  the  bibliography. 


Of  perhaps  more  general  interest  are  the  solutions  presented  by 
Reynolds  [2]   for  both  constant  and  critical  mass  flows.   These 
solutions  employ  a  number  of  simplifying  assumptions,  the  value  of  which 
are  primarily  dependent  upon  the  accuracy  of  independently  determining 
the  heat  transfer  conductances  between  the  gas,  thermal  mass,  and 
environs.  The  assumption  is  made  that  the  heat  transfer  resistances 
are  uniform  over  all  interior  and  exterior  heat  transfer  surfaces  and 
either  invariant  or  independently  predictable  with  time  depending  upon 
application.  Critical  flow  blowdown  experiments  have  been  performed 
using  time  average  free  convection  relationships  which  closely  approxi- 
mate the  conductances  predicted  from  these  solutions  [5 J  . 

Little  other  published  information  is  available  for  similar  charg- 
ing and  blowdown  experiments.   In  receivers  densely  packed  with  heat 
capacitors,  forced  convection  heat  transfer  correlations  may  apply. 
The  fluid  Jet  during  charging  may  strongly  Influence  the  convection  heat 
transfer  mechanism.   Imperfect  mixing  may  occur  during  charging  and 
cause  temperature  gradients  within  the  system. 


2.   Objectives 

Controlled  single  gas  receiver  charging  and  blowdown  experiments 
were  conducted  at  various  mass  flow  rates  and  over  varying  ranges  of 
thermal  capacitance  in  order  to  accomplish  the  following  objectives: 

(a)  Evaluate  previously  developed  simplified  analytical  expressions 
for  the  thermodynamic  state  of  the  gas  as  a  function  of  time. 

(b)  Determine  the  mechanism  of  heat  transfer  between  the  gas  in  the 
receiver  and  the  receiver  walls  during  these  processes. 

(c)  Investigate  analytical  methods  of  predicting  heat  transfer  con- 
ductances for  use  with  the  expressions  for  the  thermodynamic  state  of 
the  gas. 

(d)  Determine  the  feasibility  of  attaining  an  adiabatic  process  within 
a  reasonable  charging  or  blowdown  time. 

(e)  Investigate  the  temperature  distribution  in  the  receiver  as 
affected  by  mixing  during  charging. 


3.  Description  of  Test  Apparatus 

3.1.  General  Description 

The  test  apparatus  consisted  of  two  180  gallon  pressure  vessels  of 
150  psig  test,  one  of  which  contained  an  internal  liner  of  5/8  in. 
California  redwood.  The  tanks  were  arranged  as  shown  in  Fig.  1  with 
charging  air  introduced  at  the  bottom  through  various  size  orifice  plates. 
Each  tank  was  fitted  at  the  top  with  alternative  piping  arrangements  for 
either  critical  or  constant  mass  flow  blowdown  as  shown  in  Fig.  2. 
Physical  dimensions  are  contained  in  Table  1,  page  6   .   Air  for  charging 
the  tanks  was  supplied  from  an  air  bank  through  a  pressure  regulator 
set  at  200  psig.  Metal  strips  could  be  inserted  vertically  in  the  tanks 
to  serve  as  thermal  capacitors.   Pressure  and  temperature  were  continu- 
ously recorded  with  time.  A  photographic  view  of  the  general  system  is 
contained  in  Fig.  3. 

3.2.  Flow  Metering 

Sharp  edge  orifice  plates  were  available  in  increasing  diameters 
ranging  from  1/8  in.  through  1/2  in.  for  interchangeable  mounting  in 
flanges  attached  to  the  tank  top  and  bottom.   Discharge  coefficients  for 
orifice  sizes  through  1/4  in.  diameter  were  determined  over  a  range  of 
pressure  ratios  using  a  separate  blowdown  calibration  system.  These 
coefficients  remained  fairly  constant  up  to  the  critical  pressure  ratio. 
Discharge  coefficients  for  the  larger  orifices  were  not  determined  because 
of  limitations  in  the  calibration  system.   Instead,  experimental  data 
from  actual  charging  and  blowdown  runs  was  used  to  determine  average  flow 
rates  during  constant  flow.   Constant  flow  charging  runs  up  to  100  psig 
tank  pressure  were  attainable  from  an  air  bank  supply  of  200  psig  pressure, 
at  close  to  ambient  temperature,  upstream  from  the  orifice.   Separate 
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arrangements  permitted  either  critical  or  constant  mass  flow  blowdown. 
The  tank  top  could  be  fitted  with  a  flange  arrangement  containing  a 
seated  orifice  plate  and  a  1  in.  quick  opening  gate  valve  for  critical 
flow  runs.   Also  available  for  constant  mass  flow  blowdown  was  a  1  in. 
copper  pipe  and  flange  arrangement  equipped  with  a  hand  valve  and  bourdon 
pressure  gage  upstream  from  the  orifice  plate  for  maintaining  the  desired 
pressure  ratio. 

3.3.  Temperature  Measurements 

Each  tank  was  equipped  with  six  30  gage  copper -cons tantan  thermo- 
couples. Four  of  these  were  arranged  in  series  and  spaced  at  equal 
volumes  vertically  in  the  tank  to  give  an  average  tank  air  temperature. 
The  two  remaining  junctions  were  available  for  monitoring  wall  or  thermal 
capacitor  temperature.  Separate  thermocouples  were  available  for  moni- 
toring both  constant  flow  charging  and  constant  flow  blowdown  air  stag- 
nation temperature  upstream  from  the  orifice  plates.   Tank  air  and 
capacitor  temperatures  were  recorded  either  on  a  continuously  indicating 
Brown  or  on  a  Leeds  Northrup  potentiometer. 

3.4.  Pressure  Measurements 

Pressure  was  recorded  from  the  output  of  a  strain  gage  pressure 
pickup  of  the  diaphragm  type  [1^]   .  This  pickup  consisted  of  a  seated 
diaphragm  containing  a  center  mounted  SR4-A  X  5  strain  gage,  together 
with  two  peripheral  SR4-A5  gages  as  shown  in  Fig.  4.  This  arrangement 
was  wired  to  form  a  four  gage  external  bridge  and  provided  a  double  out- 
put with  temperature  compensation.   Calibrated  sensitivity  was  30  micro- 
inches  per  in.  per  psi.  Output  from  the  pickup  was  recorded  on  a  Baldwin 
strain  recorder.   A  bourdon  gage  was  mounted  directly  on  the  tank  for 
quick  visual  inspection  of  pressure. 
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TABLE    I 
Physical  Dimensions  of  Experimental  Apparatus 

1.  Uninsulated  Tank 

Volume:    23.9  ft3 
Internal  area:   45.5  ft 
Nominal  internal  diameter:   31  in. 
Nominal  internal  height:   54  in. 
Thermal  mass:   835  lbs 

2.  Insulated  Tank 

Volume:   17.0  ft3 

2 
Internal  area:   38.0  ft 

Nominal  internal  diameter:   29.5  in. 

Nominal  internal  height:  42  in. 

3.  Redwood  Insulation 

Thickness:   5/8  in. 

3 
Density:   25  lbs/ft 

Thermal  Conductivity:   0.06  Btu/(hr  ft  °R) 

Specific  heat:   0.6  Btu/(lb  °R) 

Diffusivity:   0.004  ft2/hr 

4.  Strip  Capacitors 

Material:   24  S  Aluminum  (Alclad) 
Mass:   10.0  lbs 
Area:   108  ft2 
Gage:   0.012  in. 


4.  Test  Procedure 

Experimental  work  consisted  of  the  following  charging  and  blow- 
down  runs : 

a.  Charging  and  blowdown  of  the  Insulated  tank  under  close  to 
adiabatic  conditions  (no  internal  heat  capacitors),  but  with  some  heat 
transfer  to  and  from  the  redwood  insulation. 

b.  Charging  and  blowdown  of  the  insulated  tank  with  added  heat 
capacitors  and  at  constant  mass  flow  rates  to  provide  a  finite  capacitance 
with  negligible  inside  heat  transfer  resistance. 

c.  Charging  and  blowdown  of  the  uninsulated  tank  at  constant  flow 
rates  such  that  the  tank  walls  provided  an  isothermal  sink  and  source 
with  a  finite  inside  heat  transfer  resistance. 

In  making  a  series  of  runs,  the  air  bank  was  first  charged  with 
290  psig  air  from  a  reciprocating  air  compressor  and  allowed  to  cool  to 
ambient  temperature.  The  air  bank  was  then  drained  of  any  condensed 
water,  and  the  pressure  regulator  set  at  200  psig.  A  quick  opening  valve 
between  the  tank  orifice  and  pressure  regulator  was  then  opened  for  a 
constant  mass  flow  charge  through  the  selected  orifice  plate.   The  piping 
system  between  the  air  bank  and  tank  was  of  sufficient  thermal  capacitance 
that  the  air  entered  at  nearly  ambient  temperature  despite  the  blowdown 
process  of  the  air  bank.  After  reaching  the  desired  tank  pressure  of 
about  100  psig,  the  tank  was  immediately  blown  down  through  either  the 
critical  flow  quick  opening  valve  or  the  constant  flow  hand  valve  arrange- 
ment.  The  latter  process  consisted  of  continually  controlling  the  hand 
valve  upstream  of  the  blowdown  orifice  plate  such  that  the  upstream 
pressure  remained  constant  at  30  psig.   Once  again,  the  piping  and  valve 
arrangement  was  of  sufficient  thermal  capacitance  that  the  temperature 
upstream  from  the  orifice  remained  nearly  ambient. 
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Throughout  a  typical  cycle  of  charge  and  blowdown,  tank  air 
temperature  was  recorded  at  5  sec.  intervals  with  the  Brown  recorder 
and  tank  air  pressure  recorded  continuously  with  the  Baldwin  recorder. 
Rune  with  added  capacitors  were  repeated  under  identical  conditions, 
recording  the  temperature  of  the  capacitors  with  the  Brown  recorder. 
Runs  varied  in  length  from  11  to  240  seconds.   Smooth  curves  were 
drawn  through  the  recorded  data  plots.   The  flow  rate  was  computed 
either  from  the  usual  critical  flow  metering  equation  using  previously 
determined  discharge  coefficients  or  from  the  experimental  mass -time 
data. 


5.  Theoretical  Considerations 
5.1.  General 

The  objective  of  analyses  of  blowdown  and  charging  systems  is  to 
obtain  relations  expressing  the  thermodynamic  state  of  the  gas  in  the 
receiver  as  a  function  of  time.  It  is  usual  to  use  temperature  as  the 
dependent  variable  with  mass  content  as  the  independent  variable.   For 
a  constant  volume  receiver,  the  pressure,  temperature,  and  mass  content 
are  related  by  the  perfect  gas  equation  of  state: 

PV  -  M(R/m)T  (1) 

Thus  having  found  temperature  as  a  function  of  mass  the  pressure  may 
also  be  found  as  a  function  of  mass,  and  knowing  the  mass -time  relation- 
ship, all  state  functions  may  be  expressed  as  functions  of  time.  We 
shall  be  concerned  in  this  investigation  primarily  with  constant  mass  flow. 

General  charging  and  blowdown  solutions  have  been  presented  by 
Reynolds  [2]  using  the  following  model: 


//////// 


*  icr 


These  analyses  assume  that  the  state  of  the  gas  is  uniform  throughout 
the  receiver  and  that  the  thermal  capacitance  of  the  receiver  walls 
and  other  internal  contents  can  be  lumped  into  a  single  capacitance, 
C  ■  Mfi  c  .  The  heat  transfer  resistance  between  the  capacitance 
and  the  environment  Is  denoted  by  R^  and  the  resistance  between  the 
capacitance  and  the  air  in  the  tank  by  R^.  The  mass  of  gas  in  the  tank 
is  M  and  the  mass  flow  rate  is  w. 
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Assuming  that  KoQ     ,  R.,  and  C  are  invariant  in  time,  a  combination 
of  energy  balance  and  heat  transfer  rate  equations  leads  to  the  general 
equation  for  charging  and  blowdown  at  constant  mass  flow  of  the  form: 


M*ff+oCS+*M*S^r  +  e  =  o   a, 


The  constants  oC   $   ">   and   $  are  functions  of  the  system  parameters 
NTU,  NTU^  ,  and  C*  .  These  parameters  are  defined  as  follows: 

_    Cr  a)* 


NTU      - 


(3) 


*Z 

1 

"o 

Rco 

Cv 

^O 

Mc 

C 

p 

u/0   cv 


NTU     =         — ! -       ^  A)oo  (4) 

c  =    &*-£*-  =   -£_  (5) 

M0  cv  M0  cv 

The  NTU  groupings  are  the  familiar  "number  of  transfer  units"  used 
in  heat  exchanger  analysis  and  represent  dimension less  conductances 
(reciprocal  of  resistance)  and  also  serve  as  a  measure  of  the  rate 
of  heat  transfer  of  the  process.  The  C*  parameter  represents  the  ratio 
of  the  capacitance  of  the  receiver  walls  and  internal  thermal  mass  to 
the  initial  capacitance  of  the  gas<> 

While  equation  (2)  may  be  solved  by  analytical,  graphical,  and  analog 

2 
computer  methods,  the  use  of  such  solutions  for  engineering  applica- 
tions is  often  difficult  or  Impracticable.   In  the  majority  of  engineering 
applications  the  magnitude  of  the  system  parameters  NTU,  NTU  oQ   ,  and 
CQ*  may  be  such  that  more  useful  closed  form  solutions  may  be  obtained. 

1  See  Appendix  I  for  the  derivation  of  the  general  equation. 

*  An  analog  computer  method  for  the  solution  of  the  general  equation 
for  blowdown  is  presented  and  evaluated  in  Appendix  III. 
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Several  of  these  solutions  have  been  developed  by  Reynolds  Q2] 
from  basic  principles,  and  it  is  with  the  application  and  evalua- 
tion of  these  solutions  for  constant  flow  that  this  thesis  is 
concerned.   Only  the  resultant  expressions  for  temperature  as  a 
function  of  mass  content  will  be  given  in  this  section  for  the 
application  of  interest.   For  convenience,  the  derivations  of  the 
general  and  of  the  simplified  solutions  are  contained  in  Appendix  I. 
5.2.   Adiabatic  Charging  and  Blowdown 

If  there  is  no  heat  transfer  to  or  from  the  gas,  the  process 
within  the  tank  is  adiabatic „   An  adiabatic  process  might  be  approxi- 
mated during  a  very  rapid  charge  or  discharge,  or  during  a  relatively 
slow  process  with  a  receiver  lined  with  an  insulator  between  the  gas 
and  wall  capacitance.  Adiabatic  behavior  corresponds  to  the  situation 
where  C  *  -  0  and  NTU  ^  »  0,  or  when  NTU  »  0.   Since  no  heat  transfer 
is  involved,  the  system  is  independent  of  mass  flow  rate.  The  follow- 
ing simple  model  and  solutions  apply: 
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-ur 


-r-r 


/  /  /    /////// 


charging 


T*   =   kT,*~ 
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5.3.   Charging  and  Blowdown  with  Negligible  Inside  Resistance 
Systems  having  a  high  inside  conductance  (negligible  inside  heat 
transfer  resistance)  combined  with  a  finite  but  low  capacitance  are 
often  found  in  service.  The  limiting  condition  of  NTU  »  oO   implies 
that  the  capacitance  temperature  is  identical  with  the  gas  tempera- 
ture. This  might  be  the  case  with  the  charging  or  blowdown  of  a  thin 
walled  pressure  vessel  where  jet  flow  and  gas  currents  create  con- 
siderable turbulence  resulting  in  high  inside  conductance  values.  The 

following  simple  model  and  solutions  apply: 

/./  /  / /////// 


*-ur„ 
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7 — 7 — 7 — 7 — 7 7      7  T* 
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blowdown 


M  +  Cfi  x 
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k-l  +  MTU^ 


*  NTO 


<£L 


One  real  value  of  these  solutions  is  that  they  supply  definite  infor- 
mation for  determining  the  effects  of  capacitance  on  the  behavior  of 
a  system. 

5.4.   Charging  and  Blowdown  with  an  Isothermal  Sink  and  Source 
In  many  systems  the  thermal  capacitance  of  the  metal  receiver 
walls  far  exceeds  the  thermal  capacitance  of  the  gas.   In  such  a  case 
the  temperature  change  of  the  walls  is  much  less  than  that  of  the  gas, 
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and  in  the  limiting  case  of  C  *  »  oO  >  heat  transfer  occurs  only 


between  an  isothermal  sink  or  source  and  the  gas.  This  might  ordinarily 
be  the  case  with  any  relatively  thick  walled  receiver  or  a  receiver 
packed  with  thermal  capacitors.  The  following  simple  model  and 
solutions  apply: 
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The  value  of  these  solutions  is  further  enhanced  in  that  they  can  be 
arranged  to  give  stepwise  changes  in  T*  when  the  parameters  NTU  and 
T^  are  varying  with  time. 
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5.5.  Isothermal  Case 
Receivers  having  a  large  thermal  capacitance  accompanied  by  high 


(13) 
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values  of  inside  heat  transfer  conductance  may  exhibit  essentially 
isothermal  behavior.  This  might  be  the  case  with  a  blowdown  wind 
tunnel  filled  with  a  metal  matrix  for  the  purpose  of  obtaining  a 
near-isothermal  blowdown.   The  effect  is  to  make  both  C0*  and  NTU 
very  large.   The  following  simple  model  applies: 


jL-A. 


'  i  '  i 


*  *  / 


/  /  / 


//////< 


/    /  / 


\TC.-T,M 


-ur„ 


-r~r 
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The  temperature  ratio  for  both  charging  and  blowdown  is  by  definition: 

T*  -  1 


(14) 


5.6.  Application 
While  it  is  not  true  that  all  engineering  systems  involving  charging 
or  blowdown  processes  may  be  treated  with  these  special  solutions, 
many  of  them  will  closely  fit  the  limiting  conditions.   To  obtain  useful 
results  with  these  constant  resistance  solutions,  one  must  use  a  suitable 
average  NTU.   In  any  real  receiver  the  wide  variations  in  flow  rates, 
temperatures,  and  pressures  will  certainly  cause  considerable  variation 
in  heat  transfer  convective  conductance.   The  remainder  of  the  report 
will  treat  in  detail  the  problem  of  obtaining  suitable  values  of  average 
or  time  varying  conductances  for  applying  these  solutions  to  typical 
systems.  Theoretical  considerations  of  heat  transfer  conductances  will 
be  discussed  where  directly  applicable  to  the  results. 


14 


6.  Summary  of  Experimental  Results 

6.1.  Form  of  Results 

Data  for  all  the  charging  and  blowdown  runs  consisted  of  tank 
air  temperature  and  pressure  measurements  as  a  function  of  time.  From 
these  data,  P*  and  T*  were  determined  from  definition  and  M*  computed 
from  the  relationship,  P*  -  M*  T*  obtained  from  equation  (1).  All 
graphical  results  are  in  the  form  of  T*  vs.  M*  since  T*  is  mora  sensi- 
tive than  P*  to  the  effects  of  heat  transfer,,  Smooth  curves  were 
drawn  through  the  plotted  data  with  the  comparable  theoretical  solution 
indicated  by  a  dashed  line.  The  adiabatic  solutions  are  also  shown 
since  any  departures  from  these  solutions  are  an  indication  of  the 
occurrence  of  heat  transfer. 

6.2.  Uninsulated  Tank 

A  series  of  constant  flow  charging  and  blowdown  runs  were  con- 
ducted with  the  uninsulated  tank  in  the  "as  found"  condition  with  no 
additional  heat  capacitors.  The  capacitance  of  the  tank  walls  was 
of  sufficient  magnitude  that  heat  transfer  essentially  occurred  to 
and  from  an  isothermal  sink  and  source.  The  following  runs  were 
conducted: 

TABLE  2 
Summary  of  Uninsulated  Tank  Runs 


Run 
No. 

Tv£e 

Orifice 
Dia.  (in) 

(lbs/hr) 

t 
(sec) 

C0* 

NTU 

1 

Charge 

1/8 

180 

200 

300 

10.5 

2 

Charge 

3/16 

380 

100 

300 

7.5 

3 

Charge 

1/4 

711 

55 

300 

5.7 

4 

Charge 

5/16 

1140 

35 

300 

4.1 

5 

Blow 

3/16 

105 

240 

41 

5.0 

6 

Blow 

1/4 

166 

140 

41 

3.5 

7 

Blow 

5/16 

257 

80 

41 

2.9 

8 

Blow 

3/8 

375 

70 

41 

2.6 

15 


The  test  results  are  presented  in  Tables  11  and  12  of  Appendix  IV 
and  are  shown  graphically  on  Figs.  5»  6,  and  7«   Solutions  based  upon 

equations  (10)  and  (11)  to  fit  the  experimental  results  are  shown  graph- 
ically for  comparison.   The  values  of  NTU  listed  in  Table  2,  page  15» 
are  those  used  in  these  equations.   The  values  of  C*  were  obtained  direct- 
ly from  equation  (5)  based  on  a  wall  thermal  capacitance  of  C0  ■  92  Btu/°R. 
6.5.    Insulated  Tank 
Several  charging  and  blowdown  runs  were  made  with  the  insulated 
tank  without  added  heat  capacitors  in  order  to  attain,  as  nearly  as 
possible,  an  adiabatic  process  and  to  gage  the  effectiveness  of  the 
insulation.   The  following  runs  were  conducted  1 

TABLE  5 
Summary  of  Insulated  Tank  Runs 


Run 
No. 

1Z21 

Orifice 
Dia.  (in) 

w0 
(lb/hr) 

t 
(seo) 

9 

Charge 

5/16 

1150 

20 

10 

Charge 

5/8 

1520 

16 

11 

Charge 

7/16 

i960 

11 

12 

Blow 

1/4 

595 

120 

15 

Blow 

5/8 

841 

4o 

14 

Blow 

7/16 

860 

55 

15 

Blow 

1/2 

1100 

26 

The  test  results  are  presented  in  Tables  1}  and  14  of  Appendix  IV  and 
are  shown  graphically  on  Figs.  8  and  9.  The  adiabatic  solutions  from 
equations  (6)  and  (7)  are  shown  for  comparison. 

6.4.    Insulated  Tank  with  Added  Heat  Capacitors 
A  series  of  constant  mass  flow  charging  and  blowdown  runs  were  made 
with  a  set  of  vertical  strip  aluminum  heat  aapaaitors  inserted  in  the 
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insulated  tank.   Capacitor  configuration  and  air  flow  rates  were 

selected  such  that  the  air  temperature  closely  followed  that  of  the 

capacitance  for  the  case  of  heat  transfer  with  negligible  inside 

resistance.   The  following  runs  were  conducted: 

TABLE  4 

Summary  of  Insulated  Tank  Runs 
with  Added  Capacitors 


Run 

No. 

Tyjae 

Orifice 
Dia.    (in) 

w 

(lb/hr) 

t 
(sec) 

* 

c 

0 

16 

Charge 

5/16 

407 

70 

10.0 

17 

Charge 

1A 

704 

40 

10.0 

18 

Charge 

5/16 

1104 

25 

10.0 

19 

Blow 

5/16 

255 

70 

1.44 

20 

Blow 

5/8 

556 

50 

1.44 

21 

Blow 

7/16 

497 

55 

1.44 

The  test  results  are  presented  in  Tables  15  and  16  of  Appendix  IV  and 
are  shown  graphically  on  Figs.  10  and  11.   Solutions  based  on  equations 
(8)  and  (9)  to  fit  the  experimental  results  are  shown  graphically  for 
comparison.   The  values  of  C0  were  obtained  directly  from  equation  (5) 
based  on  a  thermal  capacitanoe  of  Cc  -  2.2  Btu/°R. 
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7.   Discussion  of  Results 
7.1.   Uninsulated  Tank 

7.1.1.  General 

Several  charging  and  blowdown  runs,  as  listed  in  Table  2  ,  page  l£, 
were  made  with  the  uninsulated  tank  in  the  "as  found"  condition  with 
no  additional  heat  transfer  surfaces.  The  tank  was  selected  as  a 
typical  medium  pressure  air  receiver  of  the  type  utilized  for  com- 
pressed air  systems.  The  cylindrical  tank  walls  and  dished  heads 
.were  constructed  of  l/2«  in.  steel  plate  having  a  total  thermal  capa- 
citance of  Cc  -  92.0  BTU/°R. 

7.1.2.  Charging  of  Uninsulated  Tank 

The  mass  of  air  in  the  tank  prior  to  charging  was  MQ  -  1.8  lbs. 
for  all  runs,  giving  a  value  of  CQ*  ■  300  from  equation  (5).  This 
value  is  well  above  the  figure  of  C  *  «  40  recommended  by  Reynolds [2 J 
as  the  minimum  for  infinite  capacitance  behavior.   Therefore,  the 
solution  for  charging  with  heat  transfer  to  an  isothermal  sink  would 
be  expected  to  apply.  The  capacitance  temperature  during  all  charging 
runs  remained  essentially  constant  at  533  °R  giving  a  value  of 
T  *  ■  0.995.  Charging  air  entered  at  a  constant  mass  flow  rate  at 
a  temperature  of  531  °R  giving  a  value  of  Tj*  «  0.990.   Inserting 
these  constant  values  in  equation  (10)  together  with  the  best  value 
of  NTU  fitting  the  experimental  points,  resulted  in  the  solutions  of 
T*  vs.  M*  shown  graphically  on  Fig.  5.   For  values  of  M*^>  3  it  was 
found  that  equation  (10)  yields  a  constant  value  of  T*.   Inspection 
of  the  experimental  data  shows  that  for  all  runs  a  constant  value  of 
T*  was  attained  after  about  35  seconds.  The  values  of  NTU  used  in 
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equation  (10)  were  therefore  selected  to  give  a  caloulated  T*  coincid- 
ing with  the Be  constant  experimental  values. 

A  convective  heat  transfer  conductance,  h,  oan  be  found  from  the 
definition  i  NTU  =       * 

If  the  best  single  value  of  NTU  fitting  the  experimental  points  is  sub- 
stituted in  this  definition,  the  resulting  value  of  conductance,  h, 
becomes  a  mean  convective  conductance,  h,^,  implied  oonstant  over  the 
entire  charging  run.  Values  of  1^  from  NTU  are  listed  in  Table  5  below. 

TABLE  5 

Charging  Conductances 
Uninsulated  Tank 

h  max.   h  max.   h 

Tn  ave  ave 

Run  w0  NTU  (from  NTU )        (free   conv.)        (pred.Fig.12) 

No.        (lb/hr)     (Sqn.10)      (Btu/hrft2°R)     (Btu/hrft2oR)     (Btu/hrft2<>R) 

1  180  10.5  7.2  2.1  7.7 

2  580  7.5  10.7  2.5  11.0 
5  711  5.7  15.0  2.7  14.0 
4  1140  4.1  17.5  2.9  16.5 

If  the  mechanism  of  heat  transfer  between  the  air  in  the  tank  and  the 
tank  walls  is  simple   free  convection,   McAdams      [5]   recommends  as  a 
correlation  .,  .  y 


kve^f     O.I3(GrP,)3  05) 


L 

for  10      <  Gr  Pr  <10     .     Inserting  in  equation  (15)  the  maximum  value 
of  temperature  attained  during  each  run  results  in  the  maximum  values 
of  average   convective  heat  transfer   conductance,    h       ,    listed  in  Table  5. 
Since  these   conductances  are   considerably  lower  than  those   calculated 
from  NTU,    it  is  therefore  concluded   that  the  heat  transfer  mechanism  for 
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charging  is  forced  convection.  It  is  also  apparent  that  the  conduc- 
tances vary  with  time  since  the  experimental  values  of  T*  decrease 
after  attaining  an  initial  peak  value.  This  bears  out  the  observation 
of  Reynolds  C  2  1  that  in  any  real  receiver  the  wide  variations  in 
temperature  and  pressure  would  cause  considerable  variation  in  the  heat 
transfer  convective  conductance. 

The  presence  of  forced  convection  considerably  compounds  the 


problem  of  predicting  heat  transfer  conductances.  There  is  inevitably 
present  a  jet  of  fluid  discharging  through  the  orifice  into  a  region 
of  increasing  pressure.  While  no  analytical  treatment  of  this  specific 
problem  has  been  found,  there  is  considerable  information  available 
concerning  the  behavior  of  free  jets  discharging  into  a  region  of 
fluid  at  rest.  Theory  concerning  the  transport  of  momentum  in  an  iso- 
thermal, turbulent,  free  jet  of  air  discharging  at  subsonic  velocities 
into  air  at  the  same  density  might  be  expected  to  lead  to  useful  results 
in  the  present  problem. 

Therefore,  let  us  consider  the  following  to  be  a  simplified  model 
of  the  jet  and  resulting  flow  Inside  the  tank: 


Considering  the  jet  as  originating  from  a  point  at  the  orifice,  obser- 
vations show  that  such  a  jet  spreads  out  conically  with  a  width  b 
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directly  proportional  to  the  distance  x  from  the  origin  {.£}•   The 
cone  angle  has  been  found  to  be  about  25   [7]  independent  of  the  veloc- 
ity at  the  origin,  u   .  As  the  jet  advances,  the  maximum  velocity, 

OyO 

uv  ,  at  a  point  x  on  the  centerline  diminishes  proportional  to  1/x  [rJ. 
x,o 

Fresh  masses  of  fluid  are  continually  being  drawn  in  so  that  the  mass 
flow  at  succeeding  cross  sections  is  not  the  same.   The  jet  momentum, 
however,  is  constant  since  the  pressure  is  assumed  the  same  as  that  of 
the  surrounding  fluid  [&]•      These  relationships  may  be  expressed  as 
follows  I 

b     ^     X  (16) 

L  =/=>u.X/o  fT  b  -const  08) 

These  relationships  have  been  combined  and  experimentally  verified 
by  Alexander,  et  al.,  [9]  in  the  form  of  momentum  flux  ratios  expressible 
in  terms  of  the  ratio  of  orifice  diameter  to  the  distance  x,  as  follows: 


(19) 


where  K  is  a  constant  experimentally  determined  for  various  initial 
velocities.  Since  the  orifice  discharge  region  pressure  is  always  less 

than  the  oritical  pressure,  the  air  will  have  expanded  at  the  orifice 

t 
throat  only  to  the  critical  pressure.   Therefore  the  quantity  (/»u)0  0 

will  be  a  constant  independent  of  orifice  size  during  constant  flow 

charging.   Equation  (19)  may  be  placed  in  the  form, 
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a»,0   -     ^'/l     Kx        /9'A.X  <20) 

where  >o  is  now  dependent  only  on  the  tank  pressure  and  temperature, 
assumed  constant  throughout  the  tank.  From  the  familiar  critical  flow 
metering  equation  and  from  continuity,  an  average  discharge  velocity 
was  found  to  be  1050  t/sec.  At  this  velocity,  Alexander  [9]  recommends 
K»0.134,  giving  the  following  values  of  C  from  equation  (20), 

d  (in.)  C  (ft  lb)  l/2/Sec 

1/8  750 

3/16  1125 

1/4  1500 

5/16  1690 

from  which  u_    may  be  found  at  any  point  along  the  jet  axis. 

Now  consider  the  jet  as  impinging  upon  the  dished  head  of  the  tank 
and  being  deflected,  as  shown  in  the  model,  down  along  the  walls  in  a 
concentric  annulus.   If  this  be  the  case,  then  in  order  that  there  be  no 
pressure  difference  throughout  the  tank,  the  mass  flow  upwards  in  the  jet 
must  equal  the  mass  flow  down  through  the  outer  annulus  at  any  cross 
section.  From  continuity,  the  average  velocity  in  the  outer  annulus  at 
any  distance  x  in  terms  of  the  annular  and  jet  areas  Is 

IX      =   Ax,o   ~ 

A*,*.    *'°  (2D 

For  turbulent  flow  In  concentric  annuli,  McAdams   {"5  "1  recommends  for  the 
local  convective  heat  transfer  conductance, 


R*    =  -°"  (PO'W 


where  Re  is  based  on  the  equivalent  diameter, 
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D*  -  D-b„ 


(25) 


Substituting  equations  (20)  and  (21  )  in  aquation  (22)  and  using  average 
values  for  the  Prandtl  number  gives 


K^ ze  lyk 


A 


A 


ifi. 


C 


.6 


A  xA 


(22a) 


Equation  (22a)  was  used  to  calculate  the  local  convective  cond- 
uctance, h  ,  at  equally  spaced  intervals  along  the  walls  at  various 
time  intervals  for  each  charging  run  using  the  experimental  pressure- 
temperature  data  to  evaluate  the  density.   Thus  for  each  run  a  series 
or  curves  of  the  type  shown  below  were  plotted.   These  curves  show  a 


decreasing  local  conductance,  tank  top  to  bottom.   From  these  curves,  an 
average  convective  heat  transfer  conductance,  have,  was  determined  for 
each  time  increment.   Values  of  haVe  *r«  shown  graphically  on  Fig. (12). 
In  all  cases  the  predicted  h    was  found  to  have  occurred  at  a  point 

GL  V  6 

about  |  the  total  tank  height.   For  comparison,  equation  (10)  was  solved 

*      * 
for  an  NTU  corresponding  to  each  experimental  T  vs.  M  point.   Using  these 

values  of  NTU,  an  experimental  h     was  calculated  from  the  definition 

B.VG 


of  NTU  for  each  experimental  point.   These  experimental  values  of  h 
are  also  shown  on  Fig.  (12). 


ave 
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Comparison  of  the  curves  on  Fig.  (12)  shows  a  good  correlation 
between  the  predicted  h    and  that  determined  from  the  experimental 
points.  As  expected,  increases  in  pressure  and  density  predominate  at  the 
higher  flow  rates  giving  increased  convective  conductances  and  an  Increased 
NTU.  The  predicted  values  of  hflve  follow  the  trend  of  the  experimental 
values  extremely  well  although  lagging  somewhat  in  magnitude.  At  the 
lower  flow  rates  the  average  convective  conductance  remains  almost  constant 
and  this  trend  is  not  predicted  as  well  by  equation  (22a).  At  such  flow 
rates,  however,  the  process  is  essentially  isothermal  and  therefore  the 
simple  relationship  T*  a  1  gives  adequate  results.  At  the  higher  flow 
rates  where  NTU  is  changing  with  time  equation  (22a)  can  be  used  to  predict 
NTU  for  use  in  stepwise  equation  (12)  with  results  at  least  adequate  for 
design  and  certainly  better  than  would  be  predicted  from  adiabatic  behavior. 
7.1.3.  Blowdown  of  Uninsulated  Tank 

The  mass  of  air  in  the  tank  prior  to  blowdown  was  M  ■  13.2  lbs. 
for  all  runs  yielding  a  value  of  CQ*  «  41  from  equation  (5).  This  is  well 
above  the  value  of  C  *  »  6.5  recommended  by  Reynolds   [2  J  as  a  minimum 
for  infinite  capacitance  blowdown.  Therefore,  the  solution  for  blowdown 
with  heat  transfer  from  an  isothermal  source  would  be  expected  to  apply. 
The  capacitance  temperature  during  all  blowdown  runs  remained  essentially 
constant  at  535   R.   Solutions  of  T*  vs.  M*  based  on  equation  (11)  are 
shown  graphically  on  Fig.  6  and  7  for  the  best  value  of  NTU  fitting  the 
experimental  points. 

Once  again,  solving  for  the  mean  heat  transfer  conductance  over  the 
entire  run,  h^,  from  the  definition  of  NTU  yields  the  values  listed  in 
Table  6. 
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TABLE  6 


Blowdovm  Conductances 
Uninsulated  Tank 

~ 

Run 
No. 

(lb/hr) 

NTU 
(Ec.  11) 

(from  NTU) 
(Btu/hrft2oR) 

"m 
(free   conv. 
(Btu/hrft2< 

■>h 

% 
Differenoa 

5 

105 

5.0 

1.96 

1.80 

8 

6 

166 

5.5 

2.16 

1.85 

14 

7 

257 

2.9 

2.79 

1.90 

52 

8 

575 

2.6 

3.65 

1.90 

49 

It  it  now  of  interest  to  oompare  tha  test  results  with  conductances 
predicted  on  the  basis  of  an  assumption  of  a  steady  state  turbulent  free 
oonraotion  boundary  layer  over  the  entire  heating  surface.   If  average 
properties  of  air  are  inserted  in  equation  05)»  but  tha  pressure  depend- 
ence of  the  density  retained,  we  obtain  the  equation  for  an  average 
conveotive  conduotanoe, 

have  =  .036(>2AT)!'3  05.) 

where  P  is  the  pressure  in  psia  and  A  T  is  the  temperature  difference 
between  the  air  and  the  tank  walls  in  °F.  Representative  results  from 
this  equation  for  runs  5  *nd  6  are  plotted  on  Fig.  1%  In  all  oases  the 
conductance  drops  sharply  to  zero  whan  the  air  temperature  reaches  that 
of  the  tank  wall  and  then  rises  rapidly  to  remain  relatively  constant 
over  the  remainder  of  the  run.  The  mean  predicted  conductances  with  time 
from  these  ourves  are  shown  in  Table  6  above. 

In  eaoh  case,  the  predioted  mean  conductances  are  less  than  the 
measured  mean  oonduotanoes.  At  the  slower  flow  rates  the  per  oent  dif- 
ference is  reasonably  small  considering  the  inherent  uncertainty  in 
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matching  the  analytical  curves  and  the  experimental  data.   Reynolds  [5 J 
found  differences  of  the  same  order  of  magnitude  for  similar  blowdowns 
at  critical  flow.  Runs  7  and  8,  however,  show  that  this  difference 
becomes  more  pronounced  as  the  flow  rate  Increases.  Steady  state 
conditions  still  seem  to  persist  since  the  experimental  points  for 
these  runs  (Fig.  7)  follow  the  theoretical  curves  down  to  a  point  where 
constant  flow  departure  begins.  Large  flow  rates,  however,  yield  small 
values  of  NTU  and  a  closer  approach  to  adiabatic  behavior  (Fig.  7)  so 
that  the  precise  evaluation  of  conductance  might  not  be  required.   If, 
as  Reynolds  ^1  points  out,  it  is  the  relatively  slow  blowdown  process 
that  provides  the  most  difficulty  for  the  engineer  because  of  heat  transfer, 
then  for  such  cases  steady-state  heat  transfer  correlations  apparently 
provide  a  good  approximation. 

Blowdown  systems  where  the  evaluation  of  forced  convection  con- 
ductances becomes  Important  are  common  enough.  Murphy   [l]  found  that 
for  a  blowdown  wind  tunnel  with  added  heat  capacitors,  heat  transfer 
occurred  at  a  rate  of  about  eight  times  that  predicted  for  free  convec- 
tion. This  system  differed  from  the  simple  case  considered  here  in 
that  capacitors  were  installed  in  a  section  separate  from  the  air  receiver, 
through  which  the  air  passed  during  blowdown  and  charging.  Evaluation 
of  forced  convection  conductance,  therefore,  would  be  dependent  upon 
design  characteristics  and  presumably  predictable  from  available  theory. 

For  a  tank  of  dimensions  similar  to  the  one  used  in  these  experi- 
ments, it  is  probable  that  the  blowdown  heat  transfer  mechanism  would 
be  free  convection.  Thus  the  relationships  similar  to  equation  (IS) 
would  give  good  approximations  especially  if  increased  slightly  by 

26 


good  judgment  to  allow  for  some  forced  convection  effects  at  higher 

flow  rates. 

7.2.  Insulated  Tank 

7.2.1.  General 

Several  charging  and  blowdown  runs,  as  listed  in  Table  3,  page  16, 
were  conducted  with  the  insulated  tank  in  order  to  determine  how  closely 
the  adiabatic  solutions  would  apply  to  a  system  deliberately  insulated 
from  the  containing  walls.   Choice  of  an  insulation  material  presented 
a  compromise  of  several  considerations.   Ideally,  a  material  with  both 
a  low  diffusivity  (k/y£>  c  )  and  low  thermal  storage  capacity  (  p  c  ) 
is  required  in  order  to  have  the  surface  temperature  follow  that  of 
the  air  while  the  change  of  energy  storage  remains  small.   Inasmuch  as 
the  tank  was  relatively  large  in  size,  a  material  of  sufficient  struc- 
tural rigidity  was  also  needed  for  ease  of  installation  and  dimensional 
stability  under  pressure.  Foamed  polystyrene  insulation  (Dow  Chemical 
Styrofoam,  HD2)  was  first  selected  and  installed  in  the  tank.  This 
material,  although  possessing  both  a  low  diffusivity  and  low  energy 
storage  capacity,  buckled  under  a  hydrostatic  pressure  of  50  pslg  and 
came  loose  from  the  tank  walls.  This  relatively  low  pressure  also 
reduced  the  thickness  of  the  Styrofoam  about  10  percent.  Therefore, 
redwood  was  selected  as  the  most  economical  alternative  material  despite 
a  relatively  high  density  (25  lbs /ft  )  and  resulting  high  thermal 
storage  capacity, 

7.2.2.  Charging  of  Insulated  Tank 

Reynolds  [2~\   predicts  that  the  adiabatic  charging  solution  should 
be  applicable  for  values  of  NTU<^0.25  for  all  magnitudes  of  CQ*.  Prom 
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equation  (3)  it  is  seen  that  a  low  NTU  for  a  given  area  is  attained 
by  the  combination  of  a  high  initial  flow  rate  and  a  low  heat  transfer 
conductance.  Therefore,  orifice  plates  were  selected  to  give  the 
highest  flow  rate  possible  within  a  reasonable  charging  time.  The 
results  of  the  three  runs,  as  displayed  on  Fig.  8, show  considerable 
departure  from  the  adiabatlc  solution. 

Inasmuch  as  the  flow  rates  and  areas  are  fixed,  let  us  consider 
the  magnitude  of  the  mean  convective  heat  transfer  conductance,  h  , 
required  to  give  a  maximum  NTU  *  0.25.  Solving  for  h  from  the 
definition  of  NTU  yields  the  conductance  values  listed  in  Table  7. 

TABLE  7 

Charging  Conductances  -  Insulated  Tank 

max.  h 
u        /iu9u  n         (for  NTU  -  8.25) 
N0*      (lbyhr)        (Btu/hr  ft  2°R) 

9       1130  1.26 

10  1320  1.48 

11  1960  2.19 

It  was  seen,  however,  from  the  uninsulated  tank  results  of  Table  5, 
p.  19  i  that  the  conductances  are  in  reality  an  order  of  magnitude 
greater  than  the  above  maximum  values.  This  would  then  give  a 
corresponding  value  of  NTU  -  2.5.  Even  though  Fig.  8  shows  a  nearer 
approach  to  the  adiabatlc  line  for  the  fastest  run,  it  is  conceivable 
that  at  faster  charges  the  heat  transfer  conductances  would  increase 
to  such  an  extent  as  to  cause  a  departure  from  rather  than  a  con- 
tinual approach  to  the  line. 

Reference  C23  also  gives  -an   a  criterion  for  adiabatlc  charging 
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the   condition  of  CQ   <0.04  for  NTU  reasonably  low.      Even   if 

NTU  ■  2.5  is  reasonably  small,  it  will  be  shown  in  the  next  section 

that  the  C0  of  the  insulation  increases  with  time  to  values  consider- 

ably  in  excess  of  the  maximum  0Q   »  0.04.   This  low  value  of  0o  is 

unrealistic  at  best  since  its  attainment  would  require  a  capacitance 

in  this  case  of  only  C„  -  0.086  Btu/°R  (  equation   (5)  with  Mo  -  12.8  lb), 

c 

a  low  value  for  the  best  of  insulation. 

Thus  it  would  seem  that  the  criteria  of  NTU  <Cp.25,  although 
well  within  the  range  of  free  convection  heat  transfer  conductances, 
is  unattainable  with  a  tank  where  the  charging  jet  effect  raises  the 
conductance  materially.  As  a  result,  it  was  not  possible  to  attain 
a  small  enough  NTUt  while  charging  the  insulated  tank,  to  evaluate 
fully  the  limiting  conditions  for  adiabatic  charging. 
7.2.3.   Blowdown  of  Insulated  Tank 

For  the  case  of  adiabatic  blowdown,  Reynolds  [^2J  limits  the  con- 
ditions to  NTU  ^0.08  for  all  magnitudes  of  CQ.  As  in  the  case  of 
charging,  small  values  of  NTU  were  sought  by  increasing  the  flow  rate. 
Therefore,  all  blowdown  runs  were  made  at  critical  mass  flow  in  order 
to  attain  the  greatest  value  of  initial  flow  rate.   These  runs  are  sum- 
marized in  Table  5»  Page  16*  and  shown  graphically  on  Fig.  9«  As  in 
the  case  of  charging,  Fig.  9  shows  a  considerable  departure  from  adia- 
batic behavior.   If  once  again  we  solve  for  the  mean  convective  heat 
transfer  conductance,  h^,  from  the  limiting  value  of  NTU  ■  0.08,  we 
obtain  the  values  listed  in  Table  8,  page  50. 

Results  previously  obtained  from  blowdown  runs  conducted  with  the 
insulated  tank  (Table  6,  page  25)  indicate  that  conductances  slightly 
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TABLE  8 


Blowdown  Oond 

uotanoes 

Insulated 

Tank 

Run 

No. 

(lb.?hr) 

max.  hm_ 
(Btu/hr  ft^°R) 

12 

595 

0.14 

15 

841 

0.50 

14 

860 

0.51 

15 

1100 

0.29 

greater  than  those  predicted  for  free  conveotion  might  be  expected. 
Inasmuch  as  these  free  convection  conductances  are  an  order  of  magni- 
tude greater  than  the  values  of  1^  listed  in  Table  8,  the  actual  NTU 
for  blowdown  of  the  insulated  tank  is  considerably  greater  than  the 
value  of  NTU  -  0.08. 

Several  phenomena  of  interest  are  deducible  from  Fig.  9»   It  is 
noted  that  run  12,  although  at  a  slower  flow  rate  and  presumably  at  a 
higher  NTU  than  the  other  runs,  approaches  the  adiabatic  line  early  in 
the  run.   This  might  be  explained  qualitatively  by  considering  the  test 
procedure.  Blowdown  was  started  immediately  after  charging  and  thus  the 
air  was  at  a  relatively  high  temperature  while  the  insulation  remained 
essentially  at  ambient.   Therefore,  for  a  long  blowdown  there  was  time 
for  considerable  heat  transfer  to  such  an  isothermal  sink,  dropping  the 
temperature  appreciably.  At  the  end  of  the  run  conditions  were  reversed. 
The  insulation  now  acted  as  an  isothermal  source  for  heat  transfer  from 
the  lower  temperature  air.   Thus  toward  the  end  of  the  run  the  air 
temperature  dropped  less  rapidly.   (This  effect  is  also  quite  apparent 
with  blowdown  runs  5  and  6,  Fig. 6,  conducted  with  the  uninsulated  tank.) 
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As  flow  rates  increase  and  the  run  time  decreases,  this  "apparent 
adiabatic"  effect  becomes  less  pronounced  and,  in  the  case  of  run  13, 
the  temperature  drop  appears  more  nearly  linear  and  departs  further 
from  adiabatic.  As  the  flow  rate  increases  further  this  linearity  is 
more  apparent  and,  since  NTU  decreases,  the  behavior  again  approaches 
adiabatic  as  seen  by  run  14.  A  still  further  increase  in  flow  rate 
as  with  run  15  might  be  expected  to  further  approach  the  adiabatic  line, 
but  this  was  not  the  case.  Here,  it  is  probable  that  h  was  also  increas- 
ing due  to  Increased  convection,  thus  offsetting  the  increase  in  w  and 
in  effect  increasing  NTU.  Thus  it  would  appear  that  a  limiting  approach 
to  adiabatic  behavior  is  all  that  may  be  attained  within  reasonable 
blowdown  times  with  an  arrangement  such  as  this. 

Again,  Reynolds  [2]  establishes  the  condition  of  C  *  <[0.04  and 
NTU  reasonably  small  as  a  criterion  for  adiabatic  behavior.  These  con- 
ditions might  conceivably  become  more  realistic  for  blowdown  than  for 
charging  since  the  initial  mass  of  air  in  the  tank  is  considerably 
greater  than  for  charging  and  NTU  small  in  comparison.   In  this  case, 
an  M  »  8.8  lbs.  yields  a  Cc  »  0.06  Btu/°R,  a  small  value  considering 
the  capacity  of  the  redwood  insulation  for  storage  of  energy. 

7.3.  Added  Heat  Capacitors  in  Insulated  Tank 

7.3.1.  General 

Several  charging  and  blowdown  runs,  as  listed  in  Table  4,  page  17, 
were  conducted  with  added  heat  capacitors  inserted  in  the  insulated 
tank.  The  insulated  tank  was  used  in  order  to  eliminate  heat  transfer 
from  outside  the  receiver  (NTU^  -  0)  and  to  minimize  the  effects  of 
the  tank  walls  so  that  most  of  the  heating  surface  would  be  the  heat 

3l 


oapaoitors.   Sixteen  0.012  in.  thick  aluminum  strip  capacitors,  each 
12  in.  by  40  in.,  ware  installed  in  a  vertical  and  radial  position  in 
the  tank  and  spaced  evenly  in  order  to  interfere  as  little  as  possible 

with  the  air  flow  patterns  previously  discussed 0   The  surface  area  of 

o 
the  oapaoitors  was  106  ft  inoreasing  the  total  area,  inoluding  the  tank 

2 
walls,  to  146  ft  •   Total  weight  of  the  capacitors  was  10.0  lbs.  Thus 

the  primary  effect  of  the  oapaoitance  was  to  inorease  the  heat  transfer 

area  and,  therefore,  NTU. 

7.5.2«  Charging  of  Insulated  Tank  with  Added  Oapaoitors 

The  mass  of  air  in  the  tank  prior  to  charging  was  Mq  -  1 .5  lbs., 

yielding  a  value  of  0*  -10.0  from  equation  (5).   Inasmuoh  as  this  value 

of  0  is  considerably  lower  than  the  minimum  of  Oo  «  40  recommended  for 
o 

infinite  capacitance  behavior  [2],  it  is  necessary  now  to  estimate  NTU. 
Since  an  approximation  is  justifiable,  the  conductance  results  prev- 
iously obtained  with  the  uninsulated  tank  were  used.   Insertion  of  the 
minimum  experimental  values  of  average  conductance  from  Figo  12  in  the 
definition  of  NTU  yields  the  values  of  NTU  listed  in  Table  9  below. 

TABLE  9 

Charging  Conductances 
Added  Capacitors  in  Insulated  Tank 

min.  have 

Run      w0      (Fig.  12)       NTU       Q0 
No.    (lb/hr)  (Btu/hrft2oR)   (£q.  5)   

16  407        8.0        16.9     10.0 

17  714        10.5       12.6    10.0 

18  1104        12.0        9.4     10.0 

These  values  are  each  above  the  minimum  value  of  NTU>  7  recommended  by 
Reynolds  [2]  for  infinite  NTU  behavior.   Inasmuch  as  the  temperature  of 
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the  oapaoitanoe  olosely  followed  that  of  the  air,  the  solution  for 
negligible  inside  resistance  would  be  expeoted  to  apply. 
If  NTU  -  0  and  tJ  -  1,  equation  (8)  reduces  tot 
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Equation  (8a)  solved  for  oj  -10.0  is  shown  graphically  on  Fig.  10. 
Comparison  of  this  curve  with  the  experimental  mean  plot  shows  a  con- 
siderable variance.  Inasmuch  es  equation  (8a)  is  independent  of  a  heat 
transfer  rata  (i.e.  NTU^-  0),  it  is  also  independent  of  mass  flow  rate 
and  all  runs  should  follow  the  same  curve.  There  is  some  variance  here 
as  seen  in  Fig.  10. 

Since  equation  (8a)  is  sensitive  to  small  changes  in  c£,  any  in- 
crease in  receiver  oapaoitanoe,  G0,  with  time  would  considerably  affect 
T*.  This  is  especially  true  sinoe  the  initial  oapaoitanoe  of  the  air 
(Mo  oy)  is  relatively  small.  As  was  seen  from  the  insulated  tank  runs 
without  heat  capacitors  (Fig.  5)»  there  is  considerable  heat  transfer  to 
the  insulation  during  a  typical  charge.  Therefore,  a  comparison  of  the 
increase  in  insulation  Qq   alone  during  a  run  with  that  given  by  equation 
(8a)  for  the  experimental  data  should  give  an  indication  of  the  applic- 
ability of  this  solution. 

The  temperature  distribution  in  the  redwood  as  a  function  of  time 
was  calculated  by  the  Schmidt  graphieal  method  ("LOj  for  run  17»  selected 
as  representative.  The  plot  was  based  upon  a  thermal  diffusivity  of 
0.004  ft2/hr  and  an  average  surface  conductance  of  lA  Btu/(hrft2°R)  fr0m 
the  uninsulated  tank  results.  The  results  of  the  graphioal  solution  are 
shown  in  Fig.  14a.  The  oapaoitanoe  at  the  end  of  eaoh  time  interval  was 
calculated  from  the  relationship. 
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Where  Qg  is  the  heat  stored  In  redwood  as  calculated  from  the  internal 
temperature  differences  and  AT  is  the  difference  between  the  tank  air 
and  ambient  temperatures.  The  resulting  increase  in  insulation  C0*  with 
time  is  shown  as  curve  (1)  on  Fig.  14b.  Next,  values  of  C0*  for  the 
liner  plus  the  capacitors  were  calculated  from  equation  (8a)  at  each  time 
interval  using  experimental  values  of  T*  and  M*.  These  are  shown  as 
curve  (2)  on  Fig.  14b.  The  increasing  values  of  C0*  for  the  liner  alone 
were  then  subtracted  from  the  experimental  values  of  C0*  for  the  liner 
plus  the  capacitors,  resulting  in  curve  (3).  This  curve  is  essentially 
constant  at  C0*  ■  10.0  and  thus  is  in  agreement  with  the  measured  value 
for  the  added  capacitors  independent  of  the  insulation. 

Similar  treatment  of  the  other  charging  runs  would  be  expected  to 
yield  the  same  results.  Scatter  of  the  experimental  data  may  be  attri- 
buted to  variances  in  conductance,  the  smaller  flow  rates  showing  less 
effects  of  heat  transfer  to  the  insulation.  It  is  seen  then  that  equation 
(8)  gives  good  results  for  the  case  where  the  capacitance  temperature 
follows  that  of  the  gas.  For  the  more  general  application  where  external 
heat  transfer  (NTILV  0)  enters,  the  determination  of  natural  convection 
conductances  for  the  external  tank  and  ambient  air  should  present  no 
problem. 

7.3.3.  Blowdown  of  Insulated  Tank  with  Added  Capacitors 

The  mass  of  air  in  the  tank  pricr  to  blowdown  was  Mo«  9.0  lbs., 
resulting  in  CQ*  •  1.44  from  equation  (5).  This  value  is  much  lower  than 
the  infinite  capacitance  value  of  C0*»  6.5.  As  with  charging,  the 
values  of  NTU  listed  in  Table  10  have  been  estimated  from  the  results 
listed  in  Table  6,  p.  25  ,   for  the  uninsulated  tank. 
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NTU 
(Ban.   5) 

°a 

8.6 

1.44 

8.7 

1.44 

7.6 

1.44 

TABLE  10 

Blowdown  Conductances 
Added  Capacitors  In  Insulated  Tank 

Run       wo,  ha 

No_2.     (lbs/hr)     (Btu/hr  ft2oR) 

19  2551         2.9 

20  556  5.6 

21  497  4.4 
Since  the  NTU  values  are  eaoh  above  the  minimum  value  of  NTU  >  7.5 
recommended  by  Reynolds  [2]  and  since  the  oapaoitanoe  temperature  closely 
followed  that  of  the  air,  the  solutions  for  negligible  inside  resistance 
should  apply.   If  NTU^-  0  and  T1  -1,  equation  (9)  reduces  to  i 

The  solution  of  equation  (9a)  for  0*  »  1.44  is  shown  graphically 
on  Fig.  11.   The  curve  follows  the  experimental  points  quite  well  for 
the  higher  flow  rates  with  a  gradual  departure  toward  the  end  of  the 
run.  As  in  charging,  the  system  capacitance  ratio,  0o,  would  be  ex- 
peoted  to  inorease  with  time  due  to  heat  transfer  with  the  insulation. 
Such  a  change  would  be  characterized  by  lower  values  of  T  with  time. 
In  blowdown,  however,  the  initial  capacitanoe  of  the  air  (M<,  cy )  is 
much  larger  than  for  charging.   Thus  slight  increases  in  C0  with  time 
would  have  considerably  less  effect  for  blowdown  than  for  charging. 
Relatively  fast  blowdown  runs  were  made  in  order  to  minimize  the  heat 
transfer  with  the  insulation  and,  as  shown  on  Fig.  11,  the  fastest  run 
(run  21  )  follows  best  the  theoretical  curve. 

It  would  be  expected  that  the  slower  runs  would  lie  above  this 
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curve  in  a  region  of  increased  C0*  (lower  T*).  That  this  is  not  the 
case  is  shown  by  run  19.  Since  the  heat  transfer  conductances  are  much 
smaller  for  blowdovn  than  for  charging,  the  insulation  surface  tempera- 
ture does  not  follow  that  of  the  air  and  as  a  result  the  Insulation  acts 
much  as  an  Isothermal  sink  and  source.  Thus  the  slower  runs  appear 
similar  to  run  12  shown  on  Fig.  9  with  the  same  behavior  as  discussed 
previously.  All  runs,  however,  exhibit  the  effect  of  heat  transfer  from 
the  higher  temperature  insulation  to  the  lower  temperature  air. 

Despite  the  anomalies  introduced  by  the  insulation,  it  would  appear 
that  the  blowdown  solution  for  negligible  internal  resistance  yields 
reliable  predictions  for  system  behavior  where  the  capacitance  temperature 
follows  that  of  the  gas.  As  with  charging,  the  introduction  of  external 
heat  transfer  should  present  no  problem  in  the  determination  of  heat 
transfer  conductances. 
7.4.  Mixing  Considerations 

As  stated  previously,  the  derivations  of  the  solutions  evaluated 
in  this  report  are  based  upon  the  assumption  that  the  gas  temperature  is 
uniform  throughout  the  tank  at  any  given  time.  This  assumes  that  during 
charging  the  Incoming  gas  expands  and  diffuses  instantaneously  throughout 
the  gas  in  the  receiver.  Therefore,  it  was  considered  expedient  to 
perform  an  experiment  during  a  representative  charging  run  to  determine 
the  validity  of  this  assumption.  The  test  apparatus  was  so  arranged 
(Pig.  2)  that  the  four  air  temperature  thermocouples  could  be  monitored 
individually  during  a  run.  Run  4,  selected  as  representative,  was 
repeated  under  Identical  conditions,  each  time  recording  the  output  of 
a  different  thermocouple.  The  temperatures  so  obtained  revealed  no 
variance  greater  than  0.5  °F.  This  was  Indeed  the  expected  result 
considering  the  mixing  potential  of  the  Jet  flow  described  earlier. 
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This  Is,  however,  contrary  to  findings  briefly  referred  to  by 
Reynolds  [2]  to  the  effect  that  considerable  temperature  gradients 
existed  throughout  a  charging  run.  Although  the  geometry  of  the  receiver 
was  not  described,  such  results  might  be  attained  with  a  tank  where  the 
charging  jet  cone  diameter  exceeds  that  of  the  tank.  Such  an  arrange- 
ment might  conceivably  result  in  a  piston  effect  causing  temperature 
gradients  between  upper  and  lower  layers  of  gas  in  the  tank. 
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8.  Conclusions 

The  conclusions  to  be  drawn  from  this  investigation  are  summarized 
as  follows: 

t 

(a)  Excellent  quantitative  agreement  between  the  experimental  results 
and  the  simplified  analytical  solutions  of  Reynolds  [2]  have  been 
obtained  for  both  charging  and  blowdown. 

(b)  Constant  mass  charging  directly  through  an  orifice  into  the  bottom 
of  a  receiver  of  sufficient  diameter  for  the  formation  of  a  freely  ex- 
panding  fluid  jet  results  in  forced  convection  heat  transfer  conduct- 
ances between  the  air  and  the  tank  walls.  These  conductances  increase 

as  the  air  temperature  and  pressure  change  with  time  and  vary  in  magnitude 
over  the  interior  surface* 

(c)  During  relatively  slow  constant  mass  flow  blowdowns.  the  heat 
transfer  conductances  tend  to  be  constant  and  a  time  average  free  con- 
vection conductance  gives  a  satisfactory  approximation  for  a  large  cap- 
acitance system*  This  method  becomes  less  reliable  as  the  flow  rate 
increases  and  a  more  detailed  analysis,  related  to  the  geometry  of  the 
system,  would  be  required. 

(d)  Analytical  expressions  have  been  found  to  predict  forced  convection 
conductances  during  charging  that  closely  approximate  the  experimental 
values.  At  relatively  high  mass  flow  rates,  these  expressions  can  be 
used  with  available  stepwise  solutions  as  a  basis  for  preliminary  design 
of  high  capacitance  systems.  At  relatively  low  mass  flow  rates,  the 
conductances  are  of  sufficient  magnitude  to  give  essentially  isothermal 
behavior. 

(e)  Despite  the  high  mass  flow  rates  possible  with  a  relatively  large 
size  air  receiver,  it  is  unlikely  that  an  adiabatic  charging  process 
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would  ever  be  attained  within  a  reasonable  charging  time  in  the  presence 
of  a  jet  effect  forced  convection  heat  transfer  mechanism. 

(f)  Closer  approach  to  adiabatic  behavior  is  attained  during  blowdown 
than  during  charging  for  the  same  constant  flow  rate.  An  approach  limit 
is  reached  for  blowdown,  however,  since  even  though  the  flow  rate  is 
increasing  the  heat  transfer  conductance  is  also  Increasing  until  the 
effect  of  the  latter  exceeds  that  of  the  former, 

(g)  If  the  reservoir  is  of  sufficient  diameter  for  the  formation  of 
a  fully  expanding  fluid  jet  during  charging,  there  will  be  adequate 
mixing  with  no  temperature  gradient  in  the  receiver. 
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9.   Recoramenda  t  ions 

It  is  recommended  that  the  general  problem  of  charging  and  blowdown 
be  analysed  for  solution  by  analog  computer  methods.   Despite  the 
problems  encountered  with  the  simplified  computer  setup  evaluated  in 
Appendix  III,  it  is  considered  that  the  general  problem  of  charging 
and  blowdown  is  inherently  suitable  to  treatment  by  such  methods. 
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APPENDIX  I 
Derivation  of  Charging  and   Blowdown  Solutions 

1 .        Equations  of  State 

In  all   of  the   analyses   to   follow,    the  working   fluid   is   presumed   to 
be  a  perfect  gas,   with  oonstant  specifio  heats,    having  the   following 
equations  of  state:  P  V  •    ((R/m)   T 

u  -   oy   T 
X.cpT 


^  -  cp  /  °v 


((R/m) 


0      -   c 
wp       fev 


2.  General   Differential   Equation  for  Blowdown  and   Charging 
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An  energy  balanee  for  the  gas  in  the  receiver  gives: 

dt      o*         dt 

An  energy  balance  on  the  capacitanoe  gives: 

a     -    *      -     JU< 

ir*>    fa  -     di 

Rate  equations   for   the   two  heat  transfer  resistances  give: 


a./t< 


The  change  in  energy  storage  in  the  gas  due  to  mass  flow  out 
and  changes  in  specific  energy  is: 

dt      dt       dt     ^ 
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Assuming  the  specific  heat  of  the  capacitance  is  constant,  the  change  in 
energy  storage  in  the  capacitance  is: 

Combining  the  energy  balances,  rate  equations,  and  the  perfect  gas  equations 
of  state: 
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Thus   inserting  these  relationships  we  may  write: 
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These  equations  expressed  in  terms  of  dimensionless  parameters  become: 
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These  may  be  combined  into  a  single  equation  expressed  as  follows: 
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For  constant  mass  flow,  w*  »  1,  giving: 
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Similarly,  for  charging  at  constant  mass  flow 
-  c/f*   , 
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3.  Adiabatic  Charging  and  Blowdown 
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An  energy  balance  on  the  gas  in  the  receiver  gives: 

charging:   J^jM  -  d  (  M  *)  -  M  d 'u   +  U   <J  M 

blowdown:      X    dM   -    d  ( M u)  -    tfju  -h  U     d  M 
Combining  the  energy  balance  and  the  equation  of  State: 

charging:         (JlT,    -T)  d M    -    M  J? 

blowdown:    (jl-l)   T  d M     -   M  J  T 
Expressing  in  terms  of  nondimensional  variables  and  separating: 
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Integrating  from  the  initial  values  to  conditions  at  a  later  time: 
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Charging: 
Blowdown: 
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4.  Charging  and  Blowdown  with  Negligible  Inside  Resistance  at  Constant 
Ma ss^  Plow 

1X1  "-''  A  -o 


fcO 


T. 


<o 


Uc" 


'  / 

/ 


dU,  d(Mu) 
d±        dt 

Q 

>  / 

Jc>T,  rvl 


v.«V£ 


w*(-$) 


An  energy  balance  on  the  gas  and  capacitance  combined: 
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A  rate  equation  for  the  heat  transfer  resistance  R         ,  gives: 
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Since,    for  charging:     0ro  -    d  M/j-fc   and  for  Blowdown:     uT0  =  -dtA/dt 

Then,   for  charging:      g  =    *J    ^    =  or.    4J 

for  blowdown:    ^T    -    ^X    ^  ..  -  ur0     ^X 
7t  dM    dt   '  dM 

Combining  this  relationship  with  the  energy  equation,  rate  equation, 
and  the  perfect  gas  equation  of  state,  yields  the  following  differential 
equation: 
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Expressing  in  terms  of  the  diroensionless  variables  and  separating: 
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Integrating  from  the  initial  conditions  to  conditions  at  a  later  time: 
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5.  Charging  and  Blowdown  at  Constant  Mass  Flow  with  Heat  Transfer 
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To  or  From  an  Isothermal  Source 


To-.Tc 


'  I   < 


^L 


X 
x 

X 

X 

IX 


T,  M 

dt 


M 


-ft.  ur„  =  J.    dM. 


Jtuf0  =  J,(-  -T2.) 


An  energy  balance  on  the  gas  gives: 
Charging: 


Blowdown: 
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A  rate  equation  for  the  heat  transfer  resistance  R.  gives: 
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Since:      for  charging,  ^  _,    ^^1  dt       and  for  blowdown»     UTa  =  -  dM/ '  dt 


Then: 


for  charging: 
for  blowdown: 
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Combination  of  this  relationship  with  the  energy  balance,  rate  equation, 

and  perfect  gas  equation  of  state,  gives: 
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Expressing  in  terms  of  nondimensional  variables  and  separating: 
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Integrating  from  the  initial  values  to  conditions  at  a  later  time: 
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Yields  the  stepwise  equations: 
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APPENDIX  II 
Figures  1   -  14 
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Photograph  of  Test  Apparatus 
Fig.  5 
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a.  Top  View  Showing  Wiring  Details 
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APPENDIX  III 

Solution  of  the  General  Differential  Equation  for  Constant  Mass  Flow 

Blowdown  by  Analog  Computer  Methods 
1.  Derivation  of  the  Machine  Equation 

Equation  (2a) (Appendix  I)  may  be  rearranged  in  the  case  of 
NTUqq  -  0  as  follows: 

The  constants  0(  2f  £J  are  defined  in  terms  of  the  system  parameters 
as  f o 1 lows : 

oC  -  Z  -Jt   -  NTU  (25) 

7f      -  ^U  (") 

I      ^(A-l)     *m  (27) 

Equation  (2b)  is  readily  adaptable  for  solution  by  analog  computer. 
Adopting  the  standard  nomenclature  of  Korn  fill,  the  variables  of 
equation  (2b)  may  be  expressed  in  terms  of  computer  voltages  and 
functions  as  follows: 

T    =:  aTT*  (28) 

M    =  aM  M*  (29) 

_d_   _  a     A.  -  ^    jp  oo) 

In  the  above  transfer  functions,  a~  and  a  are  scale  factors  relating 
the  machine  voltages  T  and  M  to  the  variables  T*  and  M*  and  afc  is  the 
time-scale  factor  relating  differentiation  with  respect  to  the  variable 
M*  to  the  standard  operator  p  denoting  differentiation  with  respect  to 
machine  time  t. 
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Substituting  the  transfer  functions  in  equation  (2b)  results  in  the 
following  machine  equation  for  the  general  case  of  constant  mass  flow 
b lewdown : 


(2c) 


The  initial  conditions  to  be  employed  in  the  integration  of  equation  (2b) 
are  given  by  Reynolds  f~ 2 T  as  follows: 


At  M, 


1  : 


To"1 


and 


dT*\  -  X  -  1  4  NTU  -  NTU  T( 


dM 


(31) 


2.  Block  Diagram  for  Ana lag  Computer  Setup 

Machine  equation  (2c)  combined  with  the  initial  conditions  is 
represented  in  block  diagram  form  for  solution  by  d-c  analog  computer 
as  follows: 

-fT0  -T„ 


^-T 


*-   -M 


CONST. 


Figure  15 

Analog  Computer  Setup  for  Solution  of  General  Equation 
for  Constant  Flow  Slowdown 
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?•  Evaluation  of  Analog  Computer  Solutions 

The  analog  computer  blook  diagram  of  Fig.  15»  suitably  time  and 
magnitude  sealed,  was  set  up  on  a  Boeing  Model  7079  Analog  Computer  with 
the  output  to  a  Reeve  servo  plotting  board.   Itwas  originally  planned 
to  compare  the  computer  blowdown  solutions  with  those  obtained  from  the 
simplified  blowdown  equations  (Chapter  5)  best  fitting  the  experimental 
data.   The  solutions,  although  capable  of  qualitatively  reproducing  the 
general  solutions  shown  in  Fig.1  of  referenoe  [2],  were  not  sensitive 
enough  to  changes  in  the  system  parameters  to  prove  reliable. 

Solutions  for  both  the  case  of  infinite  capacitance  and  the  case  of 
negligible  internal  resistance  were  attempted  with  the  oomputer.   In  the 
infinite  oapacitanoe  case  in  which  #  and  £>  both  approach  zero  (equations 
26  and  27),  the  oomputer  plot  showed  a  barely  discernable  change  over  a 
range  of  several  integer  values  of  NTU  in  equation  (25).   The  same  held 
true  for  the  case  of  negligible  internal  resistance.   Here,  since  CQ  has 
a  finite  value,  changes  in  NTU  affect  #  and  h    as  well  as  oC  .  Once  again, 
however,  quite  large  changes  in  the  value  of  NTU  had  little  effect  on  the 
displacement  of  the  computer  solution. 

This  laok  of  sensitivity  is  in  part  due  to  the  fact  that  the  exper- 
imental data  lies,  in  all  cases,  in  a  narrow  range  of  T*  and  M  close  to 
the  origin  of  the  graphical  ploto   In  this  range,  Fig.  1  of  reference  [2] 
shows  that  there  is  little  change  in  the  trend  of  the  general  solution 
despite  large  differences  in  the  magnitude  of  the  system  parameters. 
Attempts  to  enlarge  the  scale  resulted  in  an  unreliable  performance  of  the 

plotting  equipment.   This  program,  however,  could  be  of  value  where  lower 

*      * 
values  of  T  and  M  were  of  interest,  since  in  this  region  the  general 

solution  diverges  sufficiently  to  reflect  smaller  changes  in  NTU. 


61 


APPENDIX  IV 
Experimental  Data 
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Table  11 
Charging  Run  Data-  Uninsulated  ':.^A. 


Run  1 

Orifice:  1/ 

3  in.  dia, 

t*0  =  1  ;0 

lb/hr 

C0;;-  =  300 

Tf  = 

0.990 

Tj  =  0. 

995 

t 

T 

P 

(sec) 

151 

(psia) 

T* 

zl 

__£ 

0 

537 

1U.3 

1.000 

1=00 

1.00 

5 

51*8 

17.2 

1.021 

1.16 

1.11* 

10 

551 

19.5 

1.027 

1.32 

1.29 

20 

>53 

2i;.5 

1.029 

1.65 

l.oo 

30 

552 

23.9 

1.029 

1.95 

1.90 

1*0 

552 

33.7 

1.029 

2. 23 

2.22 

5o 

552 

3J.5 

1.029 

2.59 

2.60 

6o 

552 

U2.9 

1.029 

2.°0 

2.32 

30 

552 

51.8 

1.029 

3.50 

3.U0 

100 

.52 

60.7 

1.029 

U.io 

3.93 

120 

552 

69.6 

1.029 

1*.70 

h.66 

11*0 

552 

78.5 

1.029 

5.30 

5.15 

160 

552 

S6.$ 

1.029 

5.35 

$.66 

130 

52 

9U.5 

1.029 

6.33 

6.20 

200 

552 

101.8 

1.029 

6.37 

6.63 

Run  2 


Orifice : 

3/16  in.  dia 

■. 

w0  -x.3jO  lb/hr 

C*  «  300 

Ti  = 

1    ).990 

c 

Oc903 

t 

T 

P 

* 

(sec) 

Ck) 

(psia) 

T"" 

P* 

M 

0 

537 

Uu8 

1.000 

1.00 

1.00 

10 

562 

26.7 

1.01*7 

1.30 

1.72 

20 

561 

37.9 

1.01*1* 

2.56 

2.1*5 

30 

56o 

1*9.6 

1.01*2 

3.3$ 

3.21 

Uo 

559 

60.7 

l.oia 

ii.10 

3.93 

5o 

5.9 

71.0 

1.0U1 

I*.  30 

l*.6l 

6o 

559 

30.6 

1.0U 

5.1*5 

$.16 

70 

,69 

66.S 

1.01*1 

6.00 

5.76 

30 

$$9 

95.7 

1.01*1 

6.1*7 

6.20 

•^ 

559 

99.2 

i.oia 

6.70 

6.1*3 

90 

559 

102.0 

1.  )ia 

6.90 

.  : 

9$ 

559 

loU.8 

1.01*1 

7.0  - 

6.3o 

100 

559 

106.  J 

1.01a 

7.22 

6.92 
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Table  11 

(continued) 

.-•an  3 

Orif: 

Lee:  1/ij. 

in.  dia. 

w0  - 

711  lb/hr 

cS  = 

300 

lj  =  0, 

►  990 

r*  =  0.995 

t 

T 

P 

# 

# 

(sec) 

(°R) 

(psia) 

JL 

p 

K 

0 

533 

U;.3 

1.000 

1.00 

1.00 

5 

566 

20.1; 

i.o5o 

1.33 

1.32 

10 

575 

3U.1 

l. 

2.30 

2.15 

15 

572 

U2.3 

1.061 

2.39 

2.72 

20 

%9 

52.3 

1.0 

3.53 

3.3U 

25 

563 

6l.lt 

1.056 

lul5 

3.93 

30 

563 

70.2 

1.055 

Iu71i 

k.k9 

35 

567 

73.5 

i.o51j 

5.30 

5.03 

Uo 

567 

//.I! 

1.053 

5.90 

5.60 

U5 

567 

95.2 

1.G  3 

6.14i 

6.10 

50 

567 

103.  k 

1.053 

6.9   : 

6.6U 

55 

567 

111.0 

1.053 

7.50 

7.13 

Run  k 


Orifice:  5/16 

in.  dia. 

w0 

=  llliO  lb/hr 

<- 

300 

T£  = 

0.990 

Tc  =  0.995 

t 

QT 

P 

A 

(sec) 

(°&) 

(psia) 

T""' 

r 

M 

0 

532 

lii.8 

1.000 

1.00 

1.00 

2.5 

561t 

22.2 

1.060 

i.5o 

l.lil 

5.o 

57U 

29.6 

1.079 

2.00 

1.35 

7.5 

576 

35.3 

1.0)3 

2.U2 

2.U3 

10.0 

577 

ill;.  2 

1.0JU 

2.93 

2.75 

12.5 

575 

5o.3 

1.   \1 

3.Uo 

3.H; 

i£.o 

573 

53.U 

1.077 

3.9U 

3.o3 

17.5 

572 

65.2 

1.071; 

li.liO 

U.10 

20.0 

571 

72.5 

1.072 

h.  90 

li.57 

25.0 

569 

35.3 

1.070 

5.3o 

5.1*7 

3  .0 

563 

99.2 

I.069 

6.70 

6.33 

35.0 

563 

111.3 

1.069 

7.55 

7.1U 

70 


Table  12 
IxLovdoim  Run  Data  -  Uninsulated  Tank 


Run  5 

Orifice : 

3/16  in.  di 

a.           u0 

=  105  lb/nr 

Co  =  hl- 

0 

ft 

?C 

=  0.955 

t 

T 

p 

(sec) 

(°R) 

(psia) 

y 

_£ 

H* 

0 

560 

114.8 

1.000 

1.000 

1.000 

10 

55o 

109.5 

.932 

.955 

.973 

20 

543 

l  >5.5 

.970 

.920 

.948 

ko 

533 

96^6 

.951 

.360 

.906 

60 

52U 

92.5 

.936 

.306 

.362 

30 

517 

37.3 

.921 

.760 

.323 

100 

$lh 

81.9 

.916 

.714 

.773 

120 

5o3 

77.0 

.903 

.670 

.733 

31*0 

503 

72.3 

.'99 

.630 

.702 

160 

502 

63. 2 

.o9^ 

.594 

.663 

130 

5oo 

63.2 

.393 

.550 

.617 

200 

493 

^.0 

.339 

.513 

.573 

220 

496 

5U.o 

.  J  j6 

.470 

.530 

240 

h9h 

43.2 

.  -34 

.420 

.475 

Run  6 


Orifice: 

1/i;  in.  dia, 

► 

v0  =  166  11 

D/hr 

Cq  =  Ul.o 

t|  =  0.955 

t 

T 

P 

j'. 

(■sec) 

&1 

(psia) 

T 

p"'c 

J£ 

0 

56o 

113.8 

1.000 

1.000 

1.000 

5 

553 

110.2 

.938 

.970 

.982 

10 

48 

107.0 

.980 

.940 

.958 

15 

5U3 

104.2 

.970 

.16 

.915 

20 

539 

101.0 

.964 

.930 

.924 

30 

532 

96.1 

.950 

.81l5 

.390 

ko 

526 

91.0 

.936 

.300 

.355 

5o 

513 

36.5 

. '  26 

.760 

.321 

6o 

513 

31.3 

.  15 

.720 

.7  18 

30 

503 

73.0 

.  593 

.642 

.714 

100 

496 

63.3 

.566 

.639 

120 

490 

56.4 

.375 

.495 

.566 

mo 

435 

47.7 

•  j65 

.20 

.5oo 

71 


Table  12 
(continued) 


Kun  7 


Orifice:  5/16 

i  in.  dia. 

wo  = 

:  257  lb/hr 

C*  =  la.o 

xc 

•  0.970 

t 

T 

P 

(sec 

1       £*1 

(psia) 

T* 

P* 

M 

0 

561; 

llli.3 

1.000 

1.000 

1.000 

5 

55U 

109.0 

.933 

.95o 

.966 

10 

51*7 

10U.5 

.970 

.910 

.91*0 

15 

5Uo 

100.2 

.9^ 

.371; 

.913 

20 

53l* 

96.li 

.91*3 

.31*0 

. 

25 

530 

92.0 

•9h0 

.303 

.  i  \\ 

30 

525 

33.1* 

.931 

.770 

.326 

35 

520 

31*.8 

.923 

.7liO 

.302 

ko 

516 

31.2 

.915 

.703 

.771* 

kS 

512 

77.5 

.907 

.676 

.71*6 

5o 

507 

71*.  2 

.900 

.61*6 

.717 

55 

5ol* 

70.6 

.393 

.616 

.090 

6o 

5oo 

67.2 

.837 

.536 

.660 

6$ 

1*97 

6U.0 

.331 

.^j 

.633 

70 

h9h 

6o.6 
Run  3 

.376 

.523 

.603 

Orifice:  3/3 

in.  dia. 

w0 

=  375  lb/hr 

c0  =  la.o 

S 

»  0.975 

t 

T 

P 

* 

-j<- 

# 

(sec 

)           ffiL 

(PsiaJ 

T 

P 

K 

0 

51*8 

116.1 

1.000 

1.000 

1.000 

5 

51*2 

110.3 

.939 

.950 

.961 

10 

53li 

10l*.l* 

.  75 

.900 

.923 

15 

527 

93.7 

.961 

.350 

.^ 

20 

521 

93.2 

.950 

.303 

.31*5 

25 

513 

3J.2 

.935 

.760 

.313 

30 

507 

33.3 

.921* 

.717 

.776 

ko 

1*97 

7l*.3 

.006 

.6U1 

.710 

6o 

U77 

5  .3 

•869 

.1*30 

.553 

GS 

1*75 

51.3 

M 

.1*1*2 

.511 

70 

1*71* 

1*6.8 

.361* 

.l*ol* 

.hSd 

72 


Table  13 

Charging  i'iun  Data-  Insulated  Tank 

Run  9 
Orifice:  5/16  in.  dia.    w0  =  1130  Ib/hr 
Tj  =  0.998 


t 

T 

P 

(sec) 

151 

(psia) 

X 

P" 

if 

0 

532 

Hi.3 

1.000 

1.00 

1.00 

2 

565 

23.7 

1.061 

1.60 

i.5o 

k 

532 

32.1 

1.09U 

2.17 

1.93 

6 

592 

Uo.U 

1.113 

2.73 

2.U5 

3 

593 

Ud.8 

1.126 

3.30 

2.93 

10 

60J4 

57.7 

1.135 

3.90 

3.13 

12 

607 

66.0 

1.H41 

h.k6 

3.91 

m 

610 

7U.5 

i.m6 

5.03 

1.39 

16 

612 

33.1 

l.Ut9 

5.62 

U.90 

13 

613 

91.3 

1.151 

6.20 

5.39 

20 

61U 

101.0 

1.153 

6.80 

5.90 

ittm  10 


Orifice : 

3/8  in. 

dia.           w0 

=  1320  lb/hr 

T*  =  0.5 

'93 

t 

T 

P 

(sec) 

(*R) 

(psia) 

T* 

P* 

if 

0 

523 

lli,8 

1.000 

1.00 

1.00 

1 

5U5 

13& 

1.030 

1.22 

1.13 

2 

562 

23.1 

I.060 

1.56 

l.ii7 

3 

573 

28.5 

1.033 

1.93 

1.73 

U 

533 

3U.3 

1.101 

2.32 

2.11 

5 

590 

U0.7 

1.1m 

2.75 

2.U7 

6 

5914 

hl.k 

1.123 

3.20 

2.35 

j 

601 

5U.5 

1.136 

3.63 

3.214 

10 

6o5 

60.7 

I.H1U 

U.10 

3.59 

12 

608 

73.3 

1.151 

h.9$ 

Ii.30 

m 

611 

6.3 

1.15 

5.30 

5.02 

16 

612 

96.6 

1.159 

6.53 

$.6h 

75 


Table 

13 

(continued) 

Run 

11 

Orifice : 

7/16  in. 

dia.         w0 

*  i960  lb/hr 

T*  *  l.OC 

i 

t 

T 

P 

(sec) 

m 

(psia) 

T"" 

_£ 

>-* 

h 

0 

$30 

Hi.8 

1.000 

1.00 

1.00 

1 

$3h 

17. J 

1.032 

1.20 

1.16 

2 

539 

22.3 

1.065 

1.51* 

1.1*5 

3 

561 

30.3 

1.092 

2. 

1.90 

h 

576 

35.5 

1.112 

2.67 

2.U0 

5. 

588 

1*2.6 

1.126 

3.2li. 

2.33 

6 

$<?$ 

k9.k 

1.133 

3.30 

3.3U 

7 

602 

56.3 

1.11*3 

l*.l*o 

3.31* 

3. 

608 

61.6 

1.156 

h.JO 

ii.16 

9 

612 

70.1* 

1.163 

5.53 

1.76 

10 

615 

77.6 

1.166 

6.I0 

5.2U 

11 

617 

93.9 

1.163 

6.63 

5.67 

Table  lh 

lHoTKkrwn 

Run  Data- 

■  Insulated 

Tank 

,.un  12 


Orifice: 

1/U  in.  dia. 

w0 

*  395  lb/hr 

t 

(sec) 

T 

P 

(psia) 

T 

P 

_j£ 

0 

$99 

115.7 

1.000 

1.000 

1.000 

1 

571; 

.0 

.953 

. 

.335 

20 

$$$ 

35.1 

.925 

.736 

.796 

30 

533 

7iu6 

.o99 

.61*6 

.719 

Uo 

527 

66.0 

.530 

.571 

. 

5o 

513 

53.U 

•36l» 

• 

. 

60 

,12 

51.1* 

.  53 

•Ii2»5 

.522 

70 

5  Y 

1*5.1 

•  A5 

.390 

.1*62 

3 

502 

39.3 

.339 

.31*5 

.1*12 

90 

5oo 

35.5 

.^35 

.3  7 

.■5^ 

100 

1*93 

32.U 

.  32 

•  2;0 

.337 

110 

1*97 

29.1* 

.^30 

.25h 

.306 

120 

ii96 

26.  J 

.  j29 

.232 

.2)0 

74 


Table  U4 

(continued) 

Run  13 

Orifice  s 

3/3  in.   dia. 

w0 

=    'Ail  lb/hr 

t 

T 

P 

1L 

„vt 

(sec) 

(°Li) 

(psia) 

j£ 

p" 

H 

0 

535 

111.3 

1.000 

1.000 

1.000 

2.5 

576 

loi.l* 

.935 

.910 

.921l 

567 

30.1 

•91*7 

.720 

.761 

10 

551* 

Y   .0 

.930 

.700 

.705 

15 

51*1 

69.5 

.925 

.62.; 

.676 

20 

523 

60.1 

.90U 

.51*0 

.593 

25 

518 

52.3 

.386 

.1*70 

.530 

30 

5o8 

h$.$ 

.370 

.1*09 

.1*70 

35 

502 

39.7 

.858 

.357 

.1x16 

Uo 

k96 

35.2 

.3U8 

.316 

.375 

Hun  Ik 


Orifice : 

7/16  in.  dia. 

w    = 
0 

360  lb/hr 

t 

T 

P 

,, 

a 

* 

(sec) 

en 

(psia) 

T. 

P 

M 

0 

615 

93.3 

1.000 

1.000 

1.000 

2.5 

602 

;7.3 

.976 

.338 

.9li* 

5 

537 

73.6 

.9^k 

.796 

.333 

7.5 

571* 

70.6 

.933 

.715 

.766 

10 

56k 

61i.0 

.917 

.  1*7 

.707 

12.5 

$5h 

57.2 

.900 

.57. 

.61*2 

15 

510* 

5o.7 

.^ 

.513 

.530 

17.5 

536 

1*5.7 

.372 

.1*63 

.531 

20 

523 

1*1.0 

.859 

.1*15 

.1*83 

22.5 

521 

37.1 

.8U6 

.375 

.1*1*3 

25 

511* 

3U.0 

.335 

.31*1* 

.1*12 

27.5 

5o8 

31.1 

.326 

.315 

.331 

30 

503 

29.0 

.318 

.293 

.353 

32.5 

5oo 

27.0 

.313 

.273 

.336 

35 

1*97 

2U.8 

.309 

.251 

.311 

75 


-■■able  ll; 

(continued) 

Run  15 

Or  if: 

Lee:  1/2  in.  dia. 

w0  ~ 

1100  lb/hr 

t 

T 

P 

(sec) 

(°R) 

(psia) 

-)<■ 

T 

-V. 

P 

j£ 

0 

616 

im.3 

1.000 

1.000 

1.000 

2 

6oii 

i  •k.h 

.930 

.91 

.923 

k 

592 

9U.7 

.960 

.325 

.353 

6 

581 

35.0 

.9U2 

.yUo 

.7j6 

3 

569 

75.7 

.922 

.660 

.711; 

10 

556 

67.7 

.902 

.590 

.651; 

12 

51*7 

59.7 

.837 

.520 

.586 

Ik 

537 

51.7 

.872 

•U6o 

.527 

16 

530 

li7.1 

.858 

.1W 

.JU77 

18 

522 

1;2.5 

.8I16 

.370 

•U37 

20 

515 

38.3 

.835 

•33k 

.1*00 

22 

510 

3U.8 

.828 

.303 

.367 

2k 

5o5 

31.3 

.320 

.273 

.333 

26 

h99 

28.7 

.310 

.250 

.309 

Table  15 


Charging  Rue 

1  Data-  Insulated 

Tank 

vdth  Added  Capacitors 

Run  16 

Orifice:  3/16 

in.  dia. 

wo 

=  U07  lb/hr 

g'q  =  10.0 

"1 

=  1.00 

t 

T 

P 

w 

-;;- 

-:;- 

(sec 

)         (°r0 

(psia) 

t" 

p 

II 

0 

530 

1U.8 

1.000 

1.00 

1.00 

5 

533 

22.2 

1.016 

i.5o 

1.U3 

10 

5U5 

29.6 

1.029 

2.00 

1.9U 

15 

552 

37.0 

l.OliO 

2.50 

2.1;0 

20 

556 

hk.k 

i.o5o 

3.00 

2.85 

25 

561 

51.8 

i.o59 

3.50 

3.31 

30 

561; 

59.2 

1.065 

U.00 

3.79 

^ 

567 

67.1; 

1.070 

h.tf 

U.29 

Uo 

570 

7U.6 

1.07h 

5.o5 

U.71 

16 

571 

82.2 

1.075 

$.$$ 

5.17 

^0 

572 

39.5 

l.o3o 

6.05 

5.60 

60 

575 

103.0 

1.035 

6.95 

6.U0 

70 

577 

113.6 

1.090 

7.71; 

7.10 

76 


Table  2$ 
(continued) 


Hun  17 


Orifice:  1/U 

in.  dia. 

wo  s 

70U  lb/hr 

C*  =  10.0 

- 
xl 

1.00 

t 

T 

P 

(sec) 

(°R) 

(psia) 

T* 

_£ 

Jj* 

0 

527 

ll*.8 

1.000 

1.00 

1.00 

5 

Sll 

26.7 

1.026 

1.8o 

1.75 

10 

550 

33.5 

1.0ii3 

2.60 

2*h9 

IS 

556 

51.1 

1.056 

3.16 

3.27 

20 

562 

63.6 

1.065 

U.30 

ii.Oi; 

25 

566 

77.0 

1.073 

5.20 

U.8U 

30 

570 

91.0 

1.081 

6.15 

5.68 

£ 

573 

103.6 

1.088 

7.00      , 

6.ii3 

577 

113.9 

1.095 

7.70 

7.0li 

Hun  18 


Orifice:  5/16 

in.  dia. 

wos 

llOli  lb/hr 

C*  =  10.0 

1 

1.00 

t 

T 

P 

(sec) 

ffiL 

(psia) 

T 

j£ 

JL 

0 

531 

1U.8 

1.000 

1.00 

1.00 

2.5 

538 

2luO 

1.017 

1.62 

1.60 

5 

Ski 

33.3 

1.032 

2.25 

2.13 

7.5 

tt$ 

h2.9 

1.0U5 

2.90 

2.77 

10 

560 

$2.1* 

i.o55 

3.5U 

3.35 

12.5 

%$ 

62.2 

1.061; 

U.20 

3.95 

15 

563 

71*7 

1.071 

lu35 

li.53 

17.5 

571 

31.6 

1.076 

5.52 

5.13 

20 

57U 

91.6 

1.081 

6.20 

5.73 

22.5 

576 

102.2 

1.086 

6.90 

6.36 

30 

573 

112.0 

1.090 

7.56 

6.9l* 

77 


Table  16 

KLowdo-wn  Data  -Insulated  icxk  with 
"  Added  Capacitors 


Run  19 
Orifice:  $fiJ>  in.  dia.       wc  -  255  lb/hr 


c0  -  l.Uii 

t 

T 

p 

(sec) 

iSl 

L?.3?-?! 

js! 

P* 

■if 

M 

0 

573 

ni.5 

1.000 

1.000 

1.000 

5 

567 

103.7 

.990 

.930 

.91*0 

10 

561 

98.0 

.980 

.880 

.893 

15 

^6 

91.li 

.971 

.820 

.31*5 

20 

552 

87.0 

.961* 

.730 

.810 

25 

5U8 

80.8 

.955 

.725 

.760 

30 

$hh 

73.5 

.950 

.705 

.71*2 

35 

9x2 

7U.6 

.91*6 

.670 

.706 

1*0 

539 

71.0 

.91*0 

.637 

.677 

U5 

536 

67.3 

.91*7 

.603 

.m 

5o 

53U 

63.5 

.931 

.570 

.612 

55 

531 

60.2 

.926 

.51*0 

.533 

6o 

529 

55.7 

.922 

•5oo 

•5U2 

6$ 

526 

51.3 

.913 

.1*60 

.502 

70 

521; 

1*5.3 

i 

Run  20 

.911; 

.1*12 

.1*50 

Orifice: 

3/8  in.  dia.          w0 

-  356  lb/hr 

i 

C*  -  1.1*1* 

t 

T 

P 

(sec) 

isa. 

(psia) 

j£ 

j* 

j£ 

0 

577 

112.6 

1.000 

1.000 

1.000 

5 

570 

ioU.5 

.939 

.930 

.9l*o 

10 

565 

96.6 

.979 

.360 

.879 

15 

553 

88.7 

.968 

.790 

.816 

20 

553 

82.0 

*9$9 

.730 

.762 

30 

5U3 

70.2 

.91*1 

.625 

.661; 

35 

538 

65.2 

.93U 

.580 

.622 

1*0 

535 

$9.6 

.927 

.530 

.572 

U5 

531 

52.8 

.920 

.1*70 

.512 

5o 

527 

1*6.6 

.911* 

.1*15 

.151; 

78 


Table  16 
(continued) 


iiun  21 


Orifice: 

7/16  in.  dia, 

wo 

-  U97  lb/hr 

c*  -  i.kb 

t 

J 

P 

(sec) 

£H 

(psia) 

T 

j£ 

jL 

0 

573 

lll.li 

1.000 

1.000 

1.000 

2.5 

568 

io5.h 

.993 

M 

.952 

5 

565 

100.3 

.986 

.900 

.913 

7.5 

56o 

9U.5 

.973 

.850 

.368 

10 

557 

90.1 

.972 

.810 

.333 

12.5 

tth 

85.5 

.966 

.768 

.795 

15 

55o 

81.2 

.960 

.730 

.760 

17.5 

$b$ 

76.6 

.952 

.688 

.730 

20 

5U2 

72.1i 

.9U6 

.650 

.637 

22.5 

538 

67.9 

.9hO 

.610 

.650 

25 

536 

63.8 

.935 

.573 

.613 

27.5 

533 

59.6 

.930 

.536 

.577 

30 

530 

5^.6 

.925 

•5oo 

.51il 

32.5 

527 

51.2 

.920 

•U6o 

•5oo 

35  525  1+7.2  .916  mh2k  .to 


79 
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